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Abstract 

We have studied the dynamic effect of nonlinearity on lasing in disordered 
medium. The third-order nonlinearity not only changes the frequency and 
size of lasing modes, but also modifies the laser emission intensity and laser 
pulse width. When the nonlinear response time is comparable to or longer 
than the lifetime of the lasing mode, the nonlinearity changes the laser output 
through modifying the size of the lasing mode. When the nonlinear response 
is much faster than the buildup of the lasing mode, positive nonlinearity 
always extracts more laser emission from the random medium due to the 
enhancement of single particle scattering. 
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Over the past decade, there have been considerable interest in nonlinear optical processes 
in disordered media, e.g. second harmonic generation [1-3], phase conjugation [4,5]. It has 
been shown that optical nonlinearity can modify light transport in disordered medium. [6,7] 
The coherent backscattering peak exhibits additional features in the presence of third-order 
nonlinearity. [8,9] Temporal fluctuation of scattered waves in nonlinear disordered medium 
leads to instability of speckle pattern. [10] 

Recent studies illustrate that adding gain to a disordered medium leads to lasing in the 
long-lived eigenmodes. [11-15] Despite the modes with long lifetime are preferably amplified, 
their wavefunction are not modified by the presence of gain. However, nonlinearity can 
change the eigenmodes of a disordered system. In random lasers nonlinearity is large due 
to high intensity and resonance enhancement. The laser light intensity is high, especially in 
the case of pulsed pumping. The nonlinearity at the lasing frequency, which is in resonance 
of a transition of the gain material, is resonantly enhanced. 

In this paper, we study the effect of nonlinearity on lasing in disordered medium. The 
experimental evidence of the change of lasing modes is presented first. It is followed by 
numerical simulation of a model system. The algorithm is based on the finite-difference time- 
domain (FDTD) solution to Maxwell equations coupled with the rate equations of electronic 
populations. Third-order nonlinearity with finite response time is included. Our results 
demonstrate that nonlinearity not only modifies the frequency and size of the eigenmodes 
of a disordered system, but also changes laser emission intensity and laser pulse width of 
a random laser. The nonlinear response time plays a crucial role in the dynamic nonlinear 
effect on random lasing. 
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FIG. 1. Time-integrated spectra of laser emission from a PMMA sheet with dye and micropar- 
ticles. The incident pump intensities are 2.0 12 W/m 2 (a), and 4.0 12 W/m 2 (b). The top spectrum 
is shifted vertically for clarity. 

We used two types of random media in our experiments. One is poly(methyl methacry- 
late) (PMMA) sheets containing rhodamine 640 (Rh640) perchlorate dye and titanium diox- 
ide (Ti0 2 ) microparticles, the other is ZnO nanoparticles with average diameter 100 nm. 
The dye concentration in PMMA is 10 mM. The Ti02 particles have an average diameter 
of 0.4 /xm. The density of Ti0 2 particles is ~ 1.4 x 10 12 cm" 3 . The PMMA sheets are 
optically excited by the second harmonics of a mode-locked Nd:YAG laser. The pump spot 
is ~ 100 /xm in diameter at the sample surface. The spectrum of emission from the sample 
is measured by a 0.5-meter spectrometer with a liquid nitrogen cooled CCD array detec- 
tor. When pump intensity exceeds the lasing threshold, discrete lasing modes emerge in the 
emission spectrum. As the pump intensity increases, the gain spectrum of dye molecules 
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shifts towards longer wavelength. As shown in Fig. 1, additional lasing modes appear on 
the long wavelength side, while some lasing modes on the short wavelength side disappear. 
However, individual lasing modes shift towards shorter wavelength, as marked in Fig. 1. 
Typically optical gain pulls the frequencies of lasing modes towards the peak of gain spec- 
trum (frequency pulling effect). In our experiment, the lasing modes move in the opposite 
direction as the gain spectrum. This suggests the spectral shift of lasing modes is caused 
by nonlinearity. Nonlinear Kerr effect has been observed in Rh640 solution containing TiC>2 
particles. [16] The refractive index n changes linearly with light intensity /: n = no + n2l, 
where no is linear refractive index, ri2 is the nonlinear Kerr coefficient. The nonlinearity 
is contributed mainly by dye molecules. The value of n 2 in our samples is estimated to 
be — 1CT 14 to — 1CT 15 m 2 /W. Above the lasing threshold, the excitation intensity is on the 
order of 10 11 W/m 2 . Thus the change of refractive index An ranges from — 10~ 3 to — 10~ 4 . 
This can explain the decrease of lasing wavelengths as the pump intensity increases. The 
wavelength shift is in the range 0.1 - 1.0 nm, which is on the same order of the experimental 
shift. 

In addition to time-integrated spectra, we measured time-resolved lasing spectra to track 
spectral shift of lasing modes in time. Figure 2 shows a spectral-temporal image of laser 
emission from ZnO powder. The ZnO powder is optically excited by 20 ps pulses from 
a frequency-tripled mode-locked Nd:YAG laser. The pump beam is focused to a 20 urn 
spot on the sample surface. Laser emission from the powder is dispersed by a 0.3-meter 
monochromator, then directed into a Hamamatsu streak camera. Individual lasing modes 
are separated by the monochromator, and their temporal evolutions are recorded by the 
streak camera. The lasing modes at shorter wavelength are red-shifted with time, while the 
lasing modes at longer wavelength are blue-shifted with time. The modes in the center do 
not shift in wavelength. For example, the mode labeled A shifts from 387.76 nm at t = 40 
ps to 387.96 nm at t = 130 ps. The mode labeled B shifts from 390.69 nm at t= 40 ps to 
390.58 nm at t= 105 ps. The wavelength of mode C remains nearly constant in time. 

Temporal shift of lasing frequencies must result from dynamic change of the refractive 
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index of ZnO. It is known that ZnO has a large third-order nonlinearity near its band edge. 
[17,19] The values of its third-order nonlinear coefficient range from 1CT 8 to 1CT 6 esu at 
room temperature. The nonlinear response of ZnO is also fast: its nonlinear response time 
is measured to be 160 fs near band edge. [18] 

To understand the nonlinear effect on random lasing, we performed numerical simulation 
based on the finite-difference time-domain (FDTD) solution of Maxwell equations. [20] The 
model system is one-dimensional simplification of the real experiment. The binary layers are 
made of dielectric materials with dielectric constants of e± = eo and e<i = 9e , respectively. 
The thickness of the first layer, which simulates the gain medium, is a random variable: 
a = a (l + w a a), where a = 400 nm, w a describes the strength of randomness, a is random 
value in the range [-0.5,0.5]. The thickness of the second layer, which simulates the nonlinear 
medium, is b = b (l + Wbf3), where b = 100 nm, Wb is the strength of randomness, (5 is 
random value in the range [-0.5,0.5]. For simplicity, the gain component and the nonlinear 
component are separated in our model. This separation does not change qualitatively the 
nonlinear effects that we study. The dynamics of electronic populations in the gain layers 
is described by four-level rate equations, with parameters similar to those in Ref. [21] In 
the nonlinear layers, the third-order nonlinearity is described by the nonlinear polarization 



is the third-order nonlinear coefficient. The casual response function g(t — r) = 
(1/t ) exp[— (t — t)/t ] for t > r, and g(t — r) = for t < t. r is the nonlinear response 
time. To incorporate the nonlinearity with finite response time into the FDTD algorithm, 
we introduce a new function 



[22,23] 
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The differential equation for G(x, t) can be derived as 

dG(x,t) G{x,t) | \E{x,t)\ 2 

dt T T 



(3) 
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E is related to D by accounting for both linear and nonlinear polarizations: 

E ^ = D(x,t) - P L (x,t) - P NL (x,t) 

eo^oo 

According to Eq.(4), The effective dielectric constant is 

e(x,t) = e 00 + x (1) +X i3) G(x,t). (5) 

When light frequency is far from any resonant frequencies, r — > 0, the nonlinearity becomes 
instantaneous. g(t — r) = 8(t — r), G(x,t) = \E(x,t)\ 2 . The nonlinear dielectric constant 
e = e L + x^l-El 2 ) where = + is the linear dielectric constant. 
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FIG. 2. A spectral-temporal image of laser emission from ZnO powder. 

To simulate pulsed pumping, the external pumping term in the rate equations represents 
a Gaussian pulse of width 20 ps: P r (t) = 2 x 10 8 exp[-ln2 (t - to) 2 Ao] V s > where t — 10 
ps. We first turn off the nonlinearity by setting — to simulate lasing in linear random 
medium. The random structure consists of 50 layers with w a = Wb = 0.9. The radiative 
transition of the four-level system is centered at A a = 389 nm with 7 nm spectral width. We 
choose the pumping rate so that only one mode lases. The frequency of the lasing mode does 
not exhibit any noticeable change during the lasing process. Next we include the nonlinearity. 
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We choose the value of to be 1.2 x 1CT 16 m 2 /V 2 , close to the measured value for ZnO 
at room temperature. [18] Figure 3 plots temporal evolution of lasing spectra obtained by 
Fourier transform of the electric field in 2.6 ps intervals. For negative (positive) nonlinearity 
< (x^ > 0), the lasing mode first shifts towards shorter (longer) wavelength as the 
laser intensity increases, then shifts to longer (shorter) wavelength as the laser intensity 
decreases. These results can explain qualitatively the experimental observation of lasing 
frequency shift with time. As shown in Fig. 2, the laser pulses in ZnO powder rise quickly 
then decay slowly. The rising edge of the laser pulses is too short to detect temporal shift of 
the lasing frequency. However, the falling edge of the laser pulses is long enough to observe 
noticeable shift of lasing frequency with time. In the presence of population inversion, x^ 
is positive (negative) for photon energy below (above) the band gap. For the lasing modes 
whose frequencies are below (above) ZnO band edge, since > (x^ < 0), they shift 
towards shorter (longer) wavelength in the falling edge as in Fig. 3b (3a). For the lasing 
modes at the band gap, x^ ~ 0, thus their frequencies remain nearly constant with time. 
The wavelength shift in our calculation is ~ 1-2 A, close to the experimental value. 
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FIG. 3. Temporal evolution of lasing spectra calculated for (a) = —1.2 x 1CT 16 m 2 /V 2 , (b) 

X (3) = 1.2 x 1(T 16 m 2 /V 2 . 

The spatial size of lasing modes also changes in the presence of third-order nonlinearity. 
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When only one mode lases, the size of the lasing mode can be characterized by the inverse 
partition ratio r(t) = (J \E(x,t)\ 2 dx) 2 / J \E(x,t)\ i dx. Figure 4 plots r(t) of one lasing mode 
(the same one as in Fig. 3) in both linear and nonlinear cases. For this mode, when x^ > 
(•^(3) < q) ; r(t) increases (decreases) during the lasing period. We repeat the calculation for 
lasing modes in many random configurations. However, we do not find any definite relation 
between the sign of the nonlinearity and the size change of the lasing modes. Nevertheless, 
in all simulations, we do observe a noticeable alteration in the mode size. 
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FIG. 4. Spatial size of the lasing mode r(t) calculated for x^ = (solid line), 
X (3) = -1.2 x 1(T 16 m 2 /V 2 (dotted line), and x (3) = 1.2 x 1(T 16 m 2 /V 2 (dashed line). 

Optical nonlinearity not only changes the frequency and size of lasing modes, but also 
modifies the intensity and duration of laser pulses. The change of laser output is sensitive 
to the nonlinear response time r . We calculate the laser pulses by fixing the value of 
and varying r . Figure 5 plots the total laser emission energy U(t) = (1/2) / e e\E(x, t)\ 2 dx 
for = ±1.2 x 10" 16 m 2 /V 2 . We first set r equal to the experimental value (160 fs) 
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of passive ZnO near band edge [18]. For positive (negative) x^ 3 \ the laser pulse becomes 
weaker (stronger) than that in the linear case. Next we repeat the calculation with smaller 
t since the nonlinear response time is shortened in the presence of stimulated emission. As 
shown in Fig. 5, when r = 6.5 fs, the behavior of laser pulse becomes totally different. For 
positive (negative) x^ > the laser pulse becomes stronger (weaker) and longer (shorter) than 
that in the linear case. 



2.0- 
1.5- 

E 

1.0- 

^) 

0.5- 

0.0 -L— , , 1 r , ■ 

20 25 30 35 40 45 

tips) 

FIG. 5. Laser emission energy U(t) calculated for (from top to bottom) = 1-2 x 10 -16 
m 2 /V 2 and r = 6.5 fs; x {3) = -1-2 x 10~ 16 m 2 /V 2 and r = 160 fs; x (3) = 0; x (3) = 1.2 x 10~ 16 
m 2 /V 2 and r = 160 fs; x (3) = -1-2 x 10~ 16 m 2 /V 2 and r = 6.5 fs. 

From the calculation of many random lasing modes, we conclude that the effect of non- 
linearity on laser pulse intensity and width depends on the relative magnitude of two time 
scales. One is the nonlinear response time To, the other is the lifetime r c of the lasing mode 
in the passive medium. r c can be obtained by launching a pulse at the frequency of the 
lasing mode in the passive medium and observing the temporal decay of the electric field. 
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[11] When r is comparable to or longer than r c , the change of laser pulse intensity and width 
is related to the size change of the lasing mode. If the size of the lasing mode decreases 
(increases) in the presence of nonlinearity, light confinement gets better (worse). The de- 
crease (increase) of light leakage is equivalent to an increase (decrease) of the quality factor 
of the random cavity. Hence, lasing lasts longer (shorter), and laser emission is stronger 
(weaker). When r is much shorter than r c , the change of laser output depends only on the 
sign of x^ 3 \ i- e -j positive (negative) nonlinearity always extract more (less) laser emission 
from the random medium at the same pumping rate. When the nonlinear response is much 
faster than the buildup of the lasing mode, the lasing mode cannot response fast enough 
to the nonlinear refractive index change. The phase of scattered light changes quickly due 
to rapid change of refractive index with intensity. The absence of constant phase relations 
among light waves scattered by different particles undermines the interference effect. Hence, 
the effect of single particle scattering becomes dominant over the collective effect of many 
particle scattering. For x^ > (x^ < 0), the refractive index contrast of the binary layers 
increases (decreases) as the laser intensity increases. Light scattering of a single particle 
becomes stronger (weaker). The increase (decrease) of scattering strength results in more 
(less) efficient lasing, i.e., higher (lower) laser intensity and longer (shorter) lasing period. 

In conclusion, we have studied the effect of nonlinearity on lasing in disordered medium. 
The third-order nonlinearity not only changes the frequency and size of lasing modes, but also 
modifies the laser emission intensity and laser pulse duration. When the nonlinear response 
time is comparable to or longer than the lifetime of the lasing mode, the nonlinearity changes 
the laser output through modifying the size of the lasing mode or equivalently the quality 
factor of the random cavity. When the nonlinear response is much faster than the buildup of 
the lasing mode, positive nonlinearity always extracts more laser emission from the random 
medium due to the enhancement of single particle scattering. 
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and the David and Lucille Packard Foundation. 
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